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ABSTRACT 


The  Numerical  simulations-  i.e.  the  computers  solved  the  problems  by  simulating  theoretical  models-  is  part  of 
new  technology  has  emerged  long  side  pure  theory  and  experiment  during  the  last  few  decades.  Computer  simulations  can 
bridge  the  gap  between  analysis  and  experiment.  This  paper  provides  modeling  and  simulation  of  ultrasonic  system  for 
material  characterization.  The  system  includes  the  transmitter,  transmitting  transducer,  receiver,  receiving  transducer, 
material  for  characterization.  This  paper  also  provides  the  analogy  between  transducer  and  lossy  transmission  line.  The  use 
of  PSPICE  provides  facility  for  simulation  of  complex  set  of  excitation  electronics,  the  ultrasonic  transducer,  the  material 
under  investigation,  and  the  receiving  electronics.  Initially  the  analogy  between  transducer  and  lossy  transmission  line  is 
discussed.  The  properties  for  lossy  transmission  lines  were  calculated  for  5MHz  transducer.  Thus  the  velocity  and 
attenuation  for  different  materials  is  obtained.  The  experiment  was  successfully  done  and  satisfactory  results  were 
obtained. 
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The  Numerical  simulations-  i.e.  the  use  of  computers  to  solve  problems  by  simulating  theoretical  models-  is  part 
of  new  technology  that  has  taken  place  along  side  pure  theory  and  experiment  during  the  last  few  decades. 
Numerical  simulations  help  us  to  solve  problems  that  may  be  hard  for  direct  experimental  study  or  complex  for  theoretical 
analysis  [Murawski  2002].  Computer  simulations  can  bridge  the  gap  between  experiment  and  analysis. 

These  simulations  provide  a  new  direction  in  ultrasonic  measurements,  of  great  importance  to  experiment  and 
analysis.  It  has  created  a  permanent  position  in  every  aspect  of  ultrasonic  measurements.  They  are  from  basic  research  to 
engineering  design.  The  computer  experiment  is  a  new  and  potentially  powerful  tool.  By  combining  conventional  theory, 
experiment  and  computer  simulation,  helps  to  view  unsolved  aspects  of  natural  process.  These  aspects  could  often  neither 
have  been  understood  nor  reveled  by  analysis  or  experiments  alone.  PSPICE  is  the  derivative  of  the  original  SPICE 
program  includes  several  added  values  such  as  graphical  user  interface  and  an  enhanced  plotting  [Tront  2004]. 
Its  widespread  usage  attests  to  the  applicability  of  the  program  to  large  variety  of  circuit  simulation  problems.  PSPICE 
utilizes  modified  nodal  analysis  approach.  PSPICE  can  be  used  for  non-linear  dc,  non-linear  transient  and  linear  ac  analysis 
problems.  It  can  also  be  used  to  perform  noise  analysis,  temperature  analysis,  Monte-Carlo  analysis.  One  of  the  advantages 
of  PSPICE  is  its  inclusion  of  device  models  for  active  and  passive  devices. 

PSPICE  simulators  is  used  by  engineers,  designing  a  variety  of  applications  including  RF,  Signal  processing, 
Power  supply  etc  [Yadav  et  al  1991,  Mancharkar  et  al  2005 
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The  success  of  modern  electronics  is  due  to  the  use  of  simulation  tools  which  help  to  accurately  predict  system 
behavior.  This  example  can  be  extended  to  components  such  as  piezoelectric  transducers  attached  to  the  electronics. 
The  ability  to  simulate  both  piezoelectric  transducer  and  electronics  together  renders  effective  optimizations  at  system 
level,  i.e.  minimizing  size,  cost  and  power  consumption.  In  this  paper  the  simulation  of  a  combined  electronics  and 
ultrasonic  piezoelectric  transducer  is  explored.  The  environment  used  is  compatible  with  the  electronic  simulation  tool 
PSPICE.  Improvement  and  verification  of  existing  PSPICE  models  for  ultrasonic  equipment  is  described,  and  applied  in 
the  design  of  integrated  analog  electronics  for  an  ultrasonic  measurement  system.  Ultimately  simulations  are  used  to 
predict  the  response  of  a  system  without  having  to  build  it.  From  management  point  of  view,  it  is  done  to  save  time  and 
money  as  preliminary  design  have  flaws.  It  is  not  feasible  to  have  a  scaled  model;  a  mathematical  model  should  lead  to  the 
right  direction.  In  the  case  of  ultrasound,  a  mathematical  model  makes  much  more  sense. 

DESIGN  APPROACH 
Block  Diagram  and  Description 

The  block  diagram  Ultrasonic  system  is  shown  in  figure  1.  The  high  frequency  pulse  generator  generates  sharp 
radio  frequency  pulses  of  5Mhzfrequenciy  having  pulse  width  2  to  60  (o.s.  The  repetition  rate  of  the  pulses  is  1  KHz. 
This  high  frequency  pulse  is  fed  to  the  ultrasonic  transducer  (Tx)  having  resonant  frequency  5  MHz  and  12  mm  diameter. 
The  receiver-transducer  (Rx),  which  is  at  the  other  end  of  the  sample,  converts  the  received  ultrasonic  pulse  into  an 
electrical  signal.  This  signal  is  fed  to  an  amplifier  and  then  detected. 
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Figure  1:  Block  Diagram  of  Ultrasonic  System 

Pulser-Receiver  System 

The  transducer  performs  the  conversion  of  electrical  energy  into  mechanical  energy,  and  conversely,  the 
conversion  of  mechanical  energy  into  electrical  energy.  The  ultrasonic  transmitting  transducer  transmits  the  high  frequency 
sound  waves  through  the  moving  material  without  direct  contact,  through  the  use  of  coupling  medium  i.e.  the  viscous 
liquid.  The  receiving  transducer  receives  these  ultrasound  waves  and  then  analysis  is  done  on  these  waves.  These  waves 
when  pass  through  the  material  under  test,  they  get  attenuated  and  require  finite  time  to  reach  the  receiving  transducer. 
This  time  of  flight  is  used  to  calculate  the  velocity.  This  attenuation  and  velocity  is  used  for  material  characterization. 
There  is  a  fundamental  relationship  between  the  properties  of  the  transducer  and  the  techniques  have  been  developed  for 
modeling  transducer  behavior  and  hence  predicting  the  required  system  response. 

Ultrasonic  transducer  usually  consists  of  a  piezoelectric  element  and  non-piezoelectric  layers  for  encapsulation 
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and  acoustic  matching.  Points  of  interest  are  effects  of  backing  materials,  matching  layers,  piezoelectric  materials,  layer 
thickness,  electrical  matching  and  coupling  on  transducer  characteristics  like  bandwidth,  time  response  on  different 
excitation  signals  and  ring  down  behavior. 

Wave  Theory 

The  piezoelectric  phenomena  are  modeled  using  controlled  voltage  and  current  sources  [Leach  1994].  Figure  2 
(a  and  b)  shows  the  scheme  and  controlled  source  equivalent  circuit  of  the  thickness-mode  transducer  respectively. 
Here  in  scheme  (a)  A  is  the  cross  sectional  area,  1  is  the  thickness,  f  is  the  force,  u  is  the  particle  velocity,  V  is  the  electrical 
voltage,  i  is  the  electrical  current,  h  is  the  piezoelectric  constant.  The  equivalent  circuit  consists  of  the  static  capacitance  C0 
(capacitance  between  the  electrodes),  a  transmission  line  which  represents  the  mechanical  part  of  the  piezoelectric 
transducer  and  two  controlled  sources  for  coupling  between  the  electrical  and  mechanical  part  of  the  circuit. 
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Figure  2:  Scheme  (a)  Transducer  and  (b)  Equivalent  Circuit  for  Piezoelectric  Transducer 

Suppose  an  ultrasonic  pulse  travels  through  a  medium  with  a  finite  speed  c  (m/s).  This  pulse  is  refered  as  a 
disturbance.  The  medium  reacts  to  this  disturbance.  Now  for  longitudinal  wave,  this  disturbance  is  a  compression  or 
rarefaction  of  matter.  The  medium  tries  to  displace  so  as  to  return  to  its  equilibrium  state.  The  compression  and  rarefaction 
within  the  medium  is  related  to  its  density  p  (kg/m3),  and  the  restoring  force  is  related  to  the  medium's  bulk  modulus  M 
(Pa)  [Kinder  1982].  Their  relationship  to  the  speed  of  sound  is 


Now,  an  electrical  pulse  can  travel  through  an  electrical  transmission  line.  These  pulses  are  received  at  the  other 
end  of  the  line  after  a  very  short  finite  time.  This  time  helps  to  give  the  velocity  of  these  pulses.  Analogous  to  the  acoustic 
wave,  the  electrical  pulses  are  concentrations  and  rarefaction  of  electrons  within  the  transmission  line  [Cheng,  1989]. 

A  lossy  transmission  line  can  be  built  with  distributed-parameter  network  with  the  circuit  parameters  distributed 
throughout  the  line.  One  line  segment  with  the  length  D  x  can  be  approximated  with  an  electric  circuit  as  in  the  figure  3. 
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Figure  3:  Equivalent  Circuit  of  an  Element  of  a  Transmission  Line  with  a  Length  of  A  x 

This  lossy  transmission  line  model  is  described  by  four  lumped  parameters, 

•  R  is  the  resistance  in  both  conductors  per  unit  length  in  W  /m 

•  L  is  the  inductance  in  both  conductors  per  unit  length  in  H/m 

•  G  is  the  conductance  of  the  dielectric  media  per  unit  length  in  U  Im 

•  C  is  the  capacitance  between  the  conductors  per  unit  length  in  F/m 
Where  R  and  G  is  zero  under  loss  less  conditions. 

To  derive  the  parameters,  applying  Kirchhoff  s  voltage  law  on  the  circuit  in  Figure  3: 

V(x,t)-R-Ax-i(x,t)-L-Ax-^^—V(x+Ax,t)=0  (1) 

dt 

Which  can  be  written  as: 

V(x  +  Ax,t)-V(x,t)  di(x,t)  m 

 =  R-i{x,t)  +  L   \l) 

Ax  dt 

And  then  letting  Ax  -^Owe  get: 

 ^-L-L=  R  -i(x,t)  +  L    \  (3) 

ax  at 

To  get  another  equation  relating  G  and  C,  applying  Kirchhoff  s  current  law  on  the  circuit  (figure  2)  we  get: 

i(x,t)-G-AxV(x+Ax,t)-C-AxdV(x  +  Ax,t)-i(x  +  Ax,t)  =  0  (4) 

dt 

And  letting  Ax^Oin  equation  (4)  we  get: 

_3^x10=G.v(Xf0  +  c.3V|£10  (5) 
dx  dt 

The  first  order  partial  differential  equations,  (3)  and  (5)  are  called  the  general  transmission  line  equations. 
These  equations  can  be  simplified  if  the  voltage  V(x,  t)  and  the  current  i  (x,  t)  are  time-harmonic  cosine  functions  as: 

V (x,  t)  =  real(V \x)-eiM)  (6) 

i(x,t)  =  real(I(x)-ejC0')  (7) 

Where  co  is  the  angular  frequency  Using  equation  (6)  and  equation  (7),  the  general  transmission  line  equations  (3) 
and  (4)  can  be  written  as: 

 ^  =  (R  +  J0)L)-I(x)  (8) 

ax 
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dl(x) 


dx 


■  =  {G  +  j(oC)-V\x)  (9) 


These  equations  are  called  the  time-harmonic  transmission  line  equations.  Equation  (8)  and  equation  (9)  can  be 
used  to  derive  the  propagation  constant  and  the  characteristic  impedance  of  the  line.  By  differentiating  equation  (8)  and 
equation  (9)  with  respect  to  x,  we  get: 

 -^±  =  f.y\x)  (10) 

dx 

—\L  =  f-I(x)  (11) 
dx 

Where  y  is  the  propagation  constant: 

y=a+jp  =  ^{R  +  jcoL)  ■  (G + jaQ  (12) 

The  real  part  a,  of  the  propagation  constant  is  called  the  attenuation  coefficient  in  N  p/m  and  the  imaginary  part, 
(3,  is  called  the  phase  constant  of  the  line  in  rad/m. 

The  general  solution  to  the  differential  equation  (10)  is 

V'(x)  =  A  •  e(a+jP)x  +  B  •  e(a+jP)x  (13) 

The  time  dependence  into  the  equation  (13)  can  be  obtained  by  multiplying  e'*"  to  equation  (13). 

V"(x,t)  =V\x)-eJ0X  =  A-e-ax -e^-M  +B  eax  .ei(ca^x)  (14) 

Equation  (14)  describes  two  traveling  waves;  one  traveling  in  the  positive  x  direction  with  an  amplitude  A  which 
decays  at  a  rate  a  while  the  other  with  an  amplitude  B  travels  in  the  opposite  direction  with  the  same  rate  of  decay.  The 
similar  differential  equations  govern  the  acoustical  wave  propogation.  In  the  case  of  harmonic  waves,  corresponding  to 
equations  (10)  and  (1 1),  we  have  the  lossy  linearized  acoustic  plane  wave  equations 

d2^  +  k2c-p(X,t)  =  0  (15) 
dx2 

d  u(x,t)     ,  ?  x 

 hr^  +  kf,  -U(x,t)=0  (16) 

dx2 

Where  p(x,  f)  (Pa)  is  the  pressure  and  u(x,  f)  (m/s)  is  the  particle  velocity  [Kinsleret  al  1982]  equivalent  to  y,  kcis 
the  complex  wave  number  composed  of  an  attenuation  constant  a  (N  p/m)  and  a  wave  number  k  (rad/m).  The  general 
solution  of  the  wave  equation  (15)  is, 

p(x,t)  =  A*-flV'(fflf-**)  +  BeaxeiiM+kx)  (17) 
And  is  identical  to  the  solutions  obtained  in  transmission  line's  represented  by  equation  (14).  Equation  (16)  has  a 
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solution  of  the  same  form. 

The  complex  wave  number  kc  is  given  as 

CO  1 

kc=   (18) 

C  yJl+jCOT 

Where  r  is  the  relaxation  time  [Kinsleret  al  1982] 
Using  equations  (17)  and  (18),  we  get 


CO  1     Jl  +  icotf-l 

a= — rJ   i     \  <19> 

c  ^\  +  (cot)2 


And 


c  jl  +  (C0T)2 

In  order  to  unify  the  two  theories,  an  impedance  type  analogy  is  chosen.  In  this  analogy  the  mechanical  force  is 
represented  by  voltage  and  current  represents  particle  velocity.  The  characteristic  impedance  becomes  important  at  the 
boundaries  because  at  the  boundaries  continuity  conditions  have  to  be  satisfied.  Pressure  and  normal  particle  velocity  must 
be  continuous  at  the  boundaries  as  voltage  and  current  must  be  continuous  at  connections. 

For  the  lossy  transmission  line  [Cheng,  1989]  the  characteristic  impedance  Zdis: 


R  +  jCOL 
Jel  ~^G  +  jcoC 


/.    =  .         J"  (21) 


And  for  the  lossy  acoustic  medium  [Kinsleret  al  1982]  the  characteristic  acoustic  impedance  Zflis: 

Za=  p-C-^\+jCOT  (22) 

Where  p  is  the  density  of  the  medium 

Expanding  equation  (21)  and  equation  (12),  we  obtain, 

And 

Y=-4LC{—  +  —)  +  jCo4LC  (24) 
Considering  small  but  non-negligible  losses  where  R  «coL,  G  «coC  and  cox«  1,  the  second  term  of  equation 
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L 

(4.23)  is  negligible,  leaving  the  characteristic  impedance  as  J—  Similarly  from  equation  (22),  the  low  loss  acoustical 

characteristic  impedance  can  be  approximated  as  pc.  Also,  the  wave  number  k  from  equation  (20)  become  co/c. 

To  correlate  the  two  characteristic  impedances,  impedance  type  analogy  is  used  (Figure  1  b),  here  the  force, 
(not  pressure)  is  represented  by  voltage  and  current  represents  the  particle  velocity.  Thus  the  equivalence  between  the  two 
systems  is 

Ze/=Zfl-A  (25) 

Where  A  (m2)  is  the  cross-sectional  area  of  the  acoustic  beam.  Refering  the  definition  of  the  low  loss  characteristic 
impedances  equation,  following  relationships  can  be  obtained 

L  =  Ap  (26) 

C  =  -^T  (27) 
Ape 

The  real  part  of  equation  (24)  is  the  attenuation  coefficient  a 

a  =  -  VLC  (— )  +  -  VZC  (— )  (28) 
2  L      2  C 

Drawing  a  parallel  with  the  classical  theory  of  acoustic  attenuation,  CXciassical  =  CCV  +  OCtc 

Where  is  the  coefficient  of  attenuation  due  to  viscous  losses  and  OCtc  is  the  coefficient  of  attenuation  due  to 
thermal  conduction  From  equations  (26),  (27)  and  (28),  we  can  solve  for  R  and  G  to  model  the  attenuations  such  that, 

R  =  2pcAav  (29) 

G  =  ^  (30) 
pcA 

Because  the  materials  layers  used  in  this  work  have  a  low  heat  conductance,  the  loss  due  to  thermal  conduction  is 
negligible,  and  we  let  the  conductance  G  =  0. Equations  (26),  (27),  (29)  and  (30)  are  the  final  equations  required  for 
simulation  purpose. 

Simulation  Design 
Piezoelectric  Transducer 

The  model  developed  for  piezoelectric  transducer  (Tx)  (is  shown  in  figure  4)  The  block  ^represents  the 
transmission  line.  Independent  sources  V!  and  V2  are  zero  value  sources,  which  are  used  as  ammeters  in  the  circuit.  Fj  and 
F2  are  dependent  current  sources.  The  value  of  Fj  is  given  by  F}  =h  C0  x  /  (Vj),  where  /  (Vj)  is  the  current  through  Vj. 
The  voltage  across  the  dependent  voltage  source  Ej  is  given  by  Ej=  V  (4),  where  V  (4)  is  the  voltage  at  node  4  i.e.  the 
voltage  across  Cj.  The  dependent  current  source  F2  which  charges  Cy  is  given  by  F 2  =  h  y*I(V2),  where  I(V2)  is  the  current 
through  V2  and  h  is  the  ratio  of  the  piezoelectric  stress  constant  in  the  direction  of  propagation  and  the  permittivity  with 


www.tiprc.org 


editor@tjprc.org 


88 


Apurva  Y.  Gandole,  Omkar  S.  Ghatpande  &  Tushar  R.  Jadhav 


constant  strain 

Resistor  /?y  is  included  to  prevent  node  4  from  being  a  floating  node.  From  the  mechanical  side 
(i.e.  transmission  line  Tl),  the  difference  between  the  velocity  of  each  surface  normal  to  the  propagation  path,  represented 
by  the  currents  Ui  and  u2  controls  the  current  source  Fh  The  same  transducer  is  modeled  for  receiving  transducer  (Rx). 
The  piezoelectric  material  PZT-5A  whose  material  data  was  obtained  from  [Berlincourfet  al  (2001)]  and  [Kinob  (1988)], 
was  chosen  and  given  in  table  1 


Table  1:  Physical  Properties  of  Transducers  at  25  °  C 


1 

Density  (p)  (kg/nr ) 

7750  (b) 

2 

Mechanical  Q  (Qm) 

75  (b) 

3 

Sound  velocity  (c)  (m/s) 

4350  (b) 

4 

Permitivity  with  constant  strain  (ss)  (C2/Nm2) 

7.35  xlO"9  (b) 

5 

Piezoelectric  stress  constant  (e3J)  (C/m2) 

15.8  (b) 

6 

Acoustic  Impedance  (MRayl) 

33.7  (b) 

7 

Piezoelectric  Constant  (10" 12  C/N) 

d33  =  374  (a)d15=  584  (a) 

8 

Couping  factor  (K33) 

0.66  (b) 

To  evaluate  the  model,  the  model  parameters  of  Quartz  crystal  and  PZT-5A  transducers  were  calculated  using 
equations  (26,  27,  31,  32,  33,  35)  and  given  in  table  2. 


Table  2:  Model  Parameters  of  Transducers 


SNo 

1 

Diameter  (mm) 

12.7 

2 

Cross  sectional  area  (A)  (m  ) 

0.0001267 

3 

Center  frequency  (MHz) 

5MHz 

Equivalent  Lossy  Transmission  Line  Parameters 

4 

C 

53.8nF 

5 

R 

411kQ 

6 

L 

981mH 

7 

G 

0 

8 

len 

435|im 

Electrical  Section  Parameter 

9 

Static  capacitance  Co 

2.14nF 

10 

Transmitting  constant  (h)  (N/C) 

2.15x  109 

11 

Current  source  gain  (  F2  ) 

2.15x  109 

12 

Dependant  current  source  gain  (Fi ) 

4.60 

13 

Voltage  control  voltage  source  gain  (E;) 

1 

14 

Ri 

1  KQ 

15 

c, 

IF 

The  measuring  cell  is  modeled  using  lossy  transmission  line.  To  verify  the  analogies  described  above, 
experimental  validation  for  solids  are  performed.  To  model  a  lossy  acoustic  layer  as  an  electrical  transmission  line,  the 
materials  data  i.e.  the  acoustic  beam  cross  sectional  area  A  (m2),  the  density  p  (Kg/m3),  the  speed  of  sound  (m/s)  and  the 
attenuation  constant  a  (  N  p/m)  at  the  chosen  frequency  and  temperature  is  needed  as  per  equations  26,  27,  29  and  30. 
Hence  the  materials  data  was  obtained  from  standard  references;  Lid  e(a>  (1999),  Kaye(b)ef  al  (1995) ,  Stewart(c)ef  al  (1930), 
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Greenspan(dW  al  (1959),  Carl(e),  W  east (f)  (1981),  A  say<gW  al  (1967)are  given  in  table3. 


Table  3:  Physical  Properties  of  Sample  Materials  at  25°  C 


Materials 
(Liquids) 

Density 
P 

(Kg/m3) 

Sound 
Velocity 

v  (m/s) 

Viscosity  H  (c  P) 

Acoustic 
Impedance  Z 
(M  Ray  1) 

Reference 

P 

V 

n 

Z 

Distilled  water 

1000 

1497 

0.890 

1.494 

f 

d 

e 

f 

Materials 
(Solid) 

Density 
P 

(Kg/m3) 

Sound 
Velocity 

v,  (m/s) 

Young's  Modulus 
E,  (109  N/m2) 

Acoustic 
Impedance  Z, 
(M  Ray  1) 

Reference 

P 

V 

Y 

z 

Teflon 

2140 

1390 

0.5 

2.97 

a 

g 

a 

a 

Aluminium 

2700 

5100 

69 

17.33 

a 

c 

a 

a 

Steel 

7860 

5874 

200 

46.00 

b 

c 

b 

b 

Model  parameters  of  some  sample  materials  and  intermediate  layer  are  calculated  using  the  dimensions  of 
measuring  cell,  materials  data  and  assigns  to  equivalent  lossy  transmission  line  T-intermediate  layer  and  T-sample. 
Model  parameters  of  the  some  sample  materials  are  calculated  at  25  0  C  using  equations  26,  27,  29  and  30,  which  are  given 
in  table  4. 

Table  4:  Model  Parameters  of  Some  Solid  Sample  Materials  at  25  °C 


Frequency 

5  MHz 

Model  Parameters 

C 

L 

R 

G 

SNo 

Material 

uF 

mH 

KO 

U 

1 

Distilled  water 

3.52 

126.7 

0.208 

0 

2 

Teflon 

1.88 

271 

326 

0 

3 

Aluminium 

0.071 

342 

1.74 

0 

4 

Steel 

0.029 

988 

57.3 

0 

Piezoelectric  transducer  (Ti)  Sample  Piezoelectric  transducer  (Rs) 


Figure  4:  Simulation  Circuit 

Analysis 

Time  Domain:  The  received  signals  from  simulations  in  the  time  domain.  Starting  with  solid  samples,  Figure  5, 
a  and  b  shows  the  first  20(O.S  received  by  the  5MHz  transducers  with  a  solid  sample  (Teflon)  and  Liquid  sample 
(Distilled  water)  at  25  °C 
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10us 
Time 


Al:  (3.5689u,  0.000)  A2:  (0.000,  0.000)  DIFF  (A):  (3.5689u,  0.000) 
Figure  5  (a):  Complete  Transient  Received  by  5  MHz  Transducer  at  25°C  in  Teflonat  d=0.005m 


Al:  (6.678 lu,  250.272mV)  A2:  (0.000,  0.000)  DIFF  (A):  (6.678 lu,  250.272mV) 
Figure  5  (b):  Complete  Transient  Received  by  5  MHz  Transducer  at  25°C  in  Distilled  Water  at  d=0.01m 
Attenuation  Measurement 

Experimentally,  the  attenuation  in  each  sample  is  based  on  measurement  of  pulse  height  of  received  signal  by 
varying  the  distance  between  two  transducers.  The  same  procedure  can  be  applied  to  the  simulation  by  changing  the  len 
parameter  (equivalent  to  distance  in  experimental  method)  of  T  sample  model  of  measuring  cell.  Figure  6  shows  the  probe 
screen  for  measurement  of  pulse  height  using  probe  cursor.  The  received  signal  in  distilled  water  at  25  0  C  using  5MHz 
transducer  reflects  the  maximum  pulse  height  250.272mV  with  time  of  6.678 lu  s.  Similarly  the  remaining  samples  were 
tested  and  the  results  are  compared  with  the  experimental  values,  as  shown  in  table  5  (a  &  b).  The  results  clearly  indicate 
that  the  simulation  of  an  ultrasonic  sensing  device,  including  both  electronic  circuits  and  transducer  is,  possible  using 
PSPICE.  The  major  obstacle  in  applying  this  as  a  general  method  for  ultrasound  transducer  prediction  is  the  lack  of 
appropriate  materials  data.  This  is  especially  pronounced  when  looking  for  temperature  and  frequency  dependency  of  these 
data. 


Figure  6:  Probe  Screen  for  Measurement  of  Pulse  Height  Using  Probe  Cursor 
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Table  5(a):  Comparison  of  Ultrasonic  Velocities  in  Liquid  and  Solid  Samples  at  25°C 


SNo 

Sample 

Simulated 
Velocity  (m/s) 

Literature 
Value(m/s) 

Reference 

1 

Distilled  water 

1497.431 

1497 

Weastef  al  1964 

2 

Aluminium 

5143.999 

5140 

Podesta  2002 

3 

Steel 

5870.105 

5874 

Kaye  et  al  1995 

4 

Teflon 

1398.925 

1400.31 

Yawaleefa/  1995 

5 

Copper 

4750..470 

4756 

Kaye  et  al  1995 

Table  5(b):  Comparison  of  Ultrasonic  Attenuation  in  Liquids  and  Solids  Samples  at  25°C 


SNo 

Materials 

Attenuation  at  25  °C  (2MHz)(o/f2) 
xlO  15  (  N  p  m1  Hz2) 

Reference 

Simulated 

Literature 

1 

Distilled  water 

23.100 

22 

Weast  1964 

Solid  Samples 

Attenuation  (dB/m) 

Simulated 

Literature 

2 

Aluminium 

1.199 

1.151 

Kaye  et  al  1995 

3 

Steel 

6.577 

5.600 

Kaye  et  al  1995 

4 

Teflon 

0.0048 

0.0039 

Lide  1999 

CONCLUSIONS 

Analogy  between  acoustic  media  and  transmission  lines  is  reviewed  and  an  analogous  electrical  model  using 
controlled  sources  of  an  ultrasonic  transducer  is  discussed.  A  presentation  of  Simulation  model  of  a  complete  ultrasonic 
system  is  done.  The  comparison  between  received  signal  from  the  simulation  and  that  of  an  actual  measurement  in  the  time 
domain  is  done.  This  shows  that  temperature  and  frequency  dependent  parameters  which  are  relevant  to  acoustic  wave 
propagation  can  be  modeled.  The  feasibility  can  be  demonstrated  in  an  ultrasound  transducer  setup  for  materials  property 
investigation.  Comparisons  were  made  for  attenuation  and  velocity  of  sound  for  organic  liquids  and  solids.  For  these 
materials,  the  agreement  is  good.  The  simulation  tool  provides  a  way  for  prediction  of  the  received  signal  before  anything 
is  built.  Moreover,  an  ultrasonic  simulation  package  used  helps  to  amplify  and  process  the  received  ultrasonic  signals. 
This  in  turn  provides  for  the  development  of  the  associated  electronics.  Such  a  virtual  design  and  testing  procedure  saves 
time  and  money,  and  also  provides  better  understanding  of  design  failures  and  allows  modification  of  designs  more 
efficiently  and  economically. 
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